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We present a new class of string theory solutions which are conjectured to be dual to the M = 1 conifold 
theory by Klebanov and Witten coupled to non chiral massive dynamical flavors. These are introduced, in 
the Veneziano limit, by means of suitably embedded and smeared D7-branes whose full backreaction on the 
, background is taken into account. The string solutions are used to study chromoelectric charge screening 

. effects due to the dynamical flavors in the non perturbative regime of the dual gauge theories. 
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<D ' 1 Introduction 

The nonperturbative dynamics of QCD is crucially affected by the contribution of the light quarks. One of 
the most striking effects is the maximal screening of the chromoelectric charges due to fundamental pair 
J> ■ production and chromoelectric flux breaking. On the difficult road towards constructing a string dual of 

(planar) QCD or its Af = 1 supersymmetric versions, it is thus mandatory to take care of the fundamental 
fields. These are introduced by means of Nf "flavor" D-brane sources on a closed string background 
arising from the backreaction of both the "flavor" and N c ^> 1 "color" D-branes. The flavor branes 
share with the color ones the (4d) Minkowski directions and wrap transverse non compact manifolds: 
this is how the global flavor symmetry is implemented. If the flavor branes are on top of each other 
this symmetry is non abelian. This obviously interesting setup is quite hard to deal with in practice if 
the flavor brane backreaction is not neglected: the supergravity equations of motion are a set of coupled 
partial differential equations with delta function sources. An enormous simplification can be obtained if 
the number of flavors is very large (as in the Veneziano limit 7V C , Nf — > oo with Nf /N c and g 2 N c fixed): 
by homogeneously smearing the flavor branes in the transverse space, delta function sources are replaced 
by density distributions and the relevant equations to solve are ordinary differential equations in a radial 
variable (T). The price to pay is the generic breaking of the flavor symmetry group into a product of abelian 
factors. 

If the flavors are massless, the flavor branes reach the origin of the transverse space. In the smeared 
setup the origin is a very special point: since all the Nf branes overlap there, the flavor symmetry group 
is generically enhanced. A global symmetry enhancement is usually signaled by a singularity in the string 
background: indeed all the known smeared string solutions dual to massless-flavored theories have a (good) 
singularity at the origin. If the flavors are massive the corresponding D-branes extend up to a certain 
distance from the origin and in the smeared setups there is no special point where the flavor symmetry is 
fully enhanced: the singularity due to the flavors disappears and is replaced by a region where there is no 
flavor D-brane charge at allQ The smeared density distribution, then, depends on the radial coordinate and 
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stops at the boundary of that region. As a result, the string solution inside that region is just an unfiavored 
one: this is dual to the fact that for energies smaller than the flavor mass the flavors are integrated out. 

In this paper we find new explicit string solutions accounting for the full backreaction of N c D3 and 
Nf smeared D7-branes on the conifold. Each D7-brane is holomorphically embedded in a particular way 
[5) that makes the corresponding flavor multiplets massive and the theory "non chiral" (the massless 
"chiral" setup was examined in [|4) and the massive one in |5 j). The unfiavored gauge theory is the Af = 1 
conformal SU(N C ) x SU(N C ) Klebanov-Witten (KW) quiver theory (6). Each D7-brane adds a (anti) 
fundamental flavor to a node of the gauge group. The resulting gauge theory has a UV Landau pole and a 
superpotential containing only quartic couplings and flavor mass terms. 

The new string solutions are then used to get quantitative predictions on the screening effects on the 
static quark-antiquark potential in the dual field theory at strong coupling. The analysis shows that both 
the chromoelectric force between static sources and the screening length are reduced as the number of 
flavors (resp. their mass) is increased (resp. reduced). 

2 Non chiral embeddings on the conifold 

The singular conifold is a six dimensional Calabi-Yau cone whose defining equation in C 4 can be written 
as Z\ z 2 — Z3 z± = 0. The base of the conifold is the T 1 ' 1 Sasaki-Einstein manifold whose isometry group 
is SU(2) x SU(2) x U{1). The conifold metric can be written as cIsq = dr 2 + r 2 ds^ ia , where 



and ip <G [0, 4n), (fii <E [0, 2ir), 9i £ [0, n}. The complex variables Zi can be expressed in terms of the above 
coordinates. 

The low energy dynamics of N c D3-branes stuck at the tip of the conifold is described by the Klebanov- 
Witten (KW) theory (6). This is an J\f = 1 superconformal quiver theory with gauge group SU(N C ) x 
SU(N C ) and bifundamental matter fields A, B, arranged in SU (2) x SU{2) doublets, interacting through a 
quartic superpotential Wkw = e lj e kl A.BkAjBi. The SU(2) x SU{2) x U(1)r symmetric IR fixed point 
of the theory is conjectured to have a dual description in terms of type IIB string theory on an AdS§ x T 1 ' 1 
background with N c units of five-form flux through T 1,1 and constant dilaton. 

Let us now consider a particular and well studied way of adding fundamental degrees of freedom to the 
model, i.e. by means of K-symmetric D7-branes wrapping the holomorphic 4-cycle defined by an equation 
like z\ ~ Z2 = 12 O. We call this embedding "non chiral" because, on the field theory side, a D7-brane 
adds a fundamental and an antifundamental flavor to one node of the KW quiver and the flavor mass terms 
do not break the classical flavor symmetry of the massless theory. The related superpotential is J2) 

W = Wkw + hi qi[AiBx - A 2 B 2 ]qi + h 2 h[BiA x - B 2 A 2 }q 2 + k (q iqi ) 2 + m (q iqi ) , (2) 

where we have considered a Z 2 -invariant setup where there are two stacks of D7 branes adding the same 
number of fundamental degrees of freedom to both nodes. Apart from the mass terms, the superpotential 
contains only quartic couplings^ The embedding preserves a SU (2) diagonal subgroup of the conifold 
symmetries. Moreover both the scale invariance and the U(l) isometry of the conifold are explicitly 
broken when /j ^ 0. This is precisely mapped to the explicit breaking of conformal invariance and U(1)r 
symmetry driven by the flavor masses (which are then related to the fi parameter). In the massless case 
(fj, = 0) the previous symmetries are broken at the quantum level. 

Let us now consider acting with an 50(4) ~ SU(2) x SU(2) rotation on the embedding equation 
z\ — z 2 = [i. This way we obtain the generalized embedding 

pz 1 -pz 2 + qz 3 + qz 4 = (3) 



In the case of "chiral" embeddings like z\ = /i, the flavors induce cubic couplings in W. 




(1) 
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where p, q span a unit 3-sphere 

p = cos-e H ^ > } q = sm-e ( 2 (4) 

and x G [0, 47r), € [0, 27r), G [0, 7r]. We can rewrite (01 more explicitly as 

e %P e -i{x+<i>) e i{1>-<P-L-<Pz) [r x + r2 ] = e |p, e i/3 ( (5) 

where we have defined (for reasons which will become clear in the following section) r = e p , e% Pq e l P = p, 
(f3 G [0, 2tt]) and 



Ti = sin — 
2 



To = cos - 
2 



e <(x+? a ) cos ^ sin + e «W+*>0 cos Zl sin ^ 
2 2 2 2 

sm — sm e KX+<?+vi+V2) cos — cog — 



(6) 



2 2 2 2 

We can split the complex equation (0) in two conditions 

/i = -x - + ^ - -<y9 2 + 2arg(ri +T 2 ) - 2/3 + 4?rn = , 

/a = JP\T x + T 2 \ 2 -e* p * =0, (7) 

following from the phase and the (squared) modulus, respectively. These equations imply that the phase of 
the mass terms, related to j3, drives the breaking of the shift symmetry along the ijj direction (the dual of the 
U(1)r symmetry). Moreover, since \Ti + T 2 \ 2 < 1, the D7 branes are extended along the radial direction 
only up to p = p q . The energy of a string stretched along p from the origin p = — oo to p q is a parameter 
m q which is conventionally called "constituent mass" of the flavors. When p q — > — oo, m q — > 0. 

2.1 Smearing the embeddings 

Let us now consider Nf ^> 1 branes suitably distributed within the family of embeddings (0. In particular, 
we want to restore effectively the SU(2) x SU (2) x U(l) symmetry of T 1,1 which is partially broken by 
flavor branes on top of each other. Accordingly, we place branes at different values of \i 4>i 9j w i m a 
suitably normalized symmetry preserving distribution 327r 3 pg,0, x , / 3 = Nf sm9. 

The D7 branes are sources of an F\ RR field with modified Bianchi identity dF\ = —g s Q,, where 

0= f Pe,*, x ,f>(6(fi)5(J 2 )(tfiA(tf2)d0d4>dxdp, (8) 

is the density distribution form of the smeared D7-branes. The general form of ft is dictated by the sym- 
metry, the only possibility for an exact two-form preserving SU(2) x SU{2) x U(l)^ x Z 2 being [4] 

Q = Nf ^ ( sin 9 x d9 x h dip i + sin 6 2 d9 2 A dip 2 ) — ^J^- dp A (dij) + cos 6\dtpi+ cos 6 2 d(p 2 ) , (9) 

47T 47T 

where the function Nf (p) has to follow from the smeared embedding and the dot means derivative w.r.t. 
p. Notice, from eq. (0, that VL pl p is expected to be independent on 9% , tfi and so we can just pick, say, 
9\ = 2 = in the corresponding expression deduced from (0. The integral is then easy to solve and the 
final result is inClpj, = -3N f Q(p - p q )e 3( - p "~^ which, together with ©, gives Nf(p) = 3NfQ(p - 
p q )e 3 ( Pq ~ p \ The "running effective number of flavors" thus reads 

N f (p) = , ( P < Pq ); N f (p) =N f (l- e 3 ^)) , (p > Pq ) , (10) 

where the integration constants have been fixed by requiring continuity at p q ^ As we have anticipated 
in the introduction, there is an effectively unflavored region p < p q where the D7 density distribution 



3 



In the chiral case N f (p) = N f [l- e 3 ^ -3 ^! + 3p - 3p q )] for p > p q (5). 
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vanishes. In the dual gauge theory, this corresponds to the E < m q regime where the flavors are integrated 
out. Moreover, Nf(p — > oo) = Nt and Nf(p q — > — oo) = Nf, The function Nf(p) asymptotically 
approaches the constant value Nf, it has the shape of a smoothed out Heaviside function, and it is actually 
well approximated by NfQ(p — p q ) in the small mass limit p q — > — ooQ This shape reflects the fact that at 
energies E ^> m q the gauge theory effectively behaves as the massless-flavored one. 

3 The backreacted solution 

The ansatz for the (Einstein frame) metric of the backreacted flavored KW background is of the form H 

ds 2 = h~^(p)dx 2 13 + h^(p)dsl, 

e 2 g{p) e 2 f{p) 

ds\ = e 2f ^dp 2 + — — V (d6 2 + sin 2 Oidrf) + —^—{d^ + V cosM^) 2 , (11) 
D — ' 9 — ' 

1=1,2 i=l,2 

which is a deformation of the KW metric driven by the functions f(p),g(p)- The KW solution is obtained 
when f(p) = g(p) = p. The K-symmetry equations for D7 branes on this deformed background can be 
obtained from those on the KW solution, just by substituting the KW radial coordinate r with e p . The 
ansatz for the dilaton and the RR forms reads 

4> = (f>(p) , F 5 = d h^ 1 (p) A dx" A • • • A dx 3 + Hodge dual , 
Ft = 9sN /( p > + cos fli dipx + cos 6> 2 dtp 2 ) . (12) 

The latter form explicitly solves the modified Bianchi identity dFi = —g s Q- Following a general strategy 
suggested by (9l we look for a solution of the IIB supergravity plus (smeared) D7-brane DBI+WZ action. 
As discussed in ||5] [lOl the use of the above D-brane action is well justified in the smeared setups, since 
the large transverse distance between two generic branes effectively suppresses the coupling g s Nf. Along 
the same lines of JH, the solution can be found from the supersymmetric variations of the bulk fermions 
giving rise to the following set of first order differential equations 

g = e 2 f- 2 * , / = 3 - 2e 2 /- 2 « - e <f> , 

4> = 3gJV/(p) e^, h = -277rg s N c a' 2 e~ 4 ° . (13) 

47T 

These have exactly the same form as those obtained in the (chiral) massless case by (4|, modulo the 
substitution of 7Y/ with Nf(p): this is due to the fact that the supersymmetric fermionic variations only 
contain the forms F and not the dF terms. 

The expression for the dilaton can be determined by direct integration. One gets, for p > p q 

(14) 



g s N f [ 3p - e 3 ^-3p ( e 3 P - 1 ) ] ' 

where the integration constant has been adjusted to have dilaton divergence at p = 0. This point is mapped 
to the field theory UV Landau pole. When p < p q , = e^ IR = e^ Pq ^ by continuity. 

Along the same lines, we can solve the remaining equations separately for p < p q and p > p q . Then 
we demand continuity at p = p q , regularity at p — > — oo and, as in |5], we fix the remaining integration 
constants so that h(Q) = g(0) = 0. The result is the following analytic solution. 



The step function approximation was used in a SQCD-like setup with wrapped D5-branes Q to construct the massive versions 
of the unquenched solutions found in (§). 



Copyright line will be provided by the publisher 



fdp header will be provided by the publisher 



5 



For p € (— oo, Pg] we have 



h{p) = 27irg s N c a' 2 (b + 



[2e 3p " -2-6/9,]5[l- e 3p ']-^, 



-4p 



-4p, 



4a 4 



-), b 



dp* 



(15) 



P<i 



This is just a slight deformation (driven by the constant factors) of the AdS$ xT 11 background, as can be 
realized changing variable to r = e p . In the p £ (p q , 0] region we have 



6pe 6p + 2e 3p+3p - ( e 3p - 1 ) 



(l - 6p)e 6 P + 2e 3 P+ 3 Pi ( e 3 P - 2 ) + e 6 ^ 
e s = [1 - e ^]-3 (1 - 6p) e 6 P + 2e 3p+3 P" (e 3p - 2) + e 6p ' 

/•O 

h{p) = 27Trg s N c a' 2 / e" 4 ^*)^ . 



(16) 



It can be checked that for p < the above solution has small curvature provided 1 <C Nf <C N c . This 
condition also ensures that the distance between two generic flavor branes is large (in string units). 



4 Screening from holography 

The background we have just found can be used to study how the nonperturbative gauge theory dynamics 
is affected by the dynamical flavors. The gauge theory we are talking about, the KW quiver coupled to 
massive non chiral flavors, has 2 nodes with N c colors and Nf+2N C flavors each. The two gauge couplings 
run in the same direction (the beta functions have the same positive signs) and the net result is a QED-like 
theory with a UV Landau pole at a scale Ajjv- It is interesting to ask how the expected Coulomb-like 
interaction between fundamental chromoelectric charges is affected by the sea flavors. To answer this 
question we can probe the theory with an external quark-antiquark pair Q, Q with mass Mq ^> m q , 
(Mq < Ajjv)- We can then study how the static quark-antiquark potential depends on the sea quark 
parameters m q ,Nf using well known holographic maps [llj. The QQ bound state is dual to an open 
string attached to a probe flavor-brane which extends up to pq ^> p q , (pq < 0). The string bends in the 
bulk up to a minimal radial position pq. The open string embedding is chosen as t = r, y =- cr, p = p(y) 
where y E [-L/2, L/2] is in Minkowski. The quark-anti quark distance L as well as the renormalized (by 
mass subtraction) potential V are then given by 

rPQ nv n 9 r f PQ (IF f PQ i 

where F = ^/gttgyy, G = ^gttg PP and F = F{p ). 

Using the explicit backreacted solution found above we find (see figureQ]) that the static quark-antiquark 
potential is negative and has the same Coulomb-like behavior as the unflavored conformal theory. Just as 
in the chiral setup 0, |T/(i) | decreases with Nf and increases with m q . This is the expected effect of the 
screening of the color charges due to the sea quarks: the more the theory is "unquenched" (large Nf, small 
m q ) the more the modulus of the quark-antiquark force is reduced. Finally, the large L behavior of V(L) 
is not strongly affected by the choice of the cutoff Mq; in the small L region, instead, we see that 
is an increasing function of Mq. 

The above analysis refers to the connected part of the potential, i.e. that which is obtained by artificially 
turning off the mixing of the QQ state with Qq, qQ pairs which can originate from the decay of the initial 
meson. The string setup we have constructed is well aware that this decay can happen as the dynamical 
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Fig. 1 From left to right, plots of V(L) at different values of: g s Nf = 1, 0.6, 0.2, top to bottom (at m q — 0.1, 
Mq = 3); m q = 0.5,1,1.5, top to bottom (at g a N f = 1, Mq = 2.5); Mq = 0.6,1.2,3, top to bottom (at 
g s Nf = 1, m q — 0.1). We have put {2na')~ 1 = 1 and ^j27irg s N c a[' 2 — 1 in the expressions for V and L 
respectively. 

q, q pair is popped out of the vacuum. Indeed, due to the presence of the dynamical flavor branes, the open 
string attached to the D7 probe can break, causing the decay. The reason why the connected part of the 
potential is an interesting object in this and related setups is the fact that, due to the smearing over the large 
transverse volume, each decay of this kind is suppressed as l/N c and not as Nf /N c as it would happen if 
the flavor branes would be on top of each other. Hence the QQ state is metastable. 

In the present setup the smallest heavy-light mesons which can be produced are nearly massless (due to 
the smearing, there will always exist a dynamical flavor brane intersecting the probe Q one). The critical 
length for the decay to happen can be thus estimated as the solution of 2Mq + V(L S ) ~ 0. It is possible to 
verify that, just as in the chiral case Q, the screening length L s decreases with Nf and Mq and increases 
with m q . 
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